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Abstract 

We interpret lattice data for the equation of state of pure gauge 
SU{3)c by an evaporation model. At low temperatures gluons are 
frozen inside the gluon condensate, whose dynamics is described in 
terms of a dilaton lagrangian. Above the critical temperature quasi- 
free gluons evaporate from the condensate: a first order transition is 
obtained by minimizing the thermodynamical potential of the system. 
Within the model it is possible to reproduce lattice QCD results at 
finite temperature for thermodynamical quantities such as pressure 
and energy. The gluonic longitudinal mass can also be evaluated; it 
vanishes below the critical temperature, where it shows a discontinuity. 
At very large temperatures we recover the perturbative scenario and 
gluons are the only asymptotic degrees of freedom. 
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In recent years, precise data have become available concerning QCD ther¬ 
modynamics at high temperature via numerical simulation on a lattice (for 
a recent review, see P). Data exist now both for the pure gauge sector and 
for complete QCD at zero chemical potential; the latter has been explored 
both in the limit of inhnite quark masses and in the chiral limit. Moreover, 
calculations for hnite values of quark masses and for a non-vanishing chem¬ 
ical potential are now appearing; the availability of such a large number of 
new lattice data surely represents an important opportunity to test the ef¬ 
fectiveness of models in reproducing the hnite temperature phase transition. 

In this paper, we will concentrate on the simplest case, i.e., the pure 
gauge sector. The main known characteristics of SU{3)c at hnite temper¬ 
ature are the following. A hrst order transition takes place at a tempera¬ 
ture Tc = (271 ± 2) MeV and at T ~ few the asymptotic limit of a 
Stefan-Boltzmann gas is not yet reached. A small but not negligible value 
for pressure, entropy and energy at T just below Tc has been computed and 
the size of the discontinuity of the energy at Tc, representing the latent heat 
of the transition, has also been estimated. 

In our work, we interpret lattice data by assuming a theoretical scenario 
similar to the one suggested by Simonov [2]. In this approach, at T < Tc 
the dynamics of the gluon condensate is dominant, while at T > Tc the 
condensate evaporates in the form of quasi-free gluons. 

Several models have been used to describe lattice data (for a review, see 
j2j). Early attempts were based on the MIT bag model |3], but more sophis¬ 
ticated approaches became necessary when more precise lattice data started 
appearing: for instance, the results for the energy density, the pressure and 
the entropy of a pure gluon system at T > Tc can be well reproduced by quasi¬ 
particle models in which gluons acquire an effective temperature-dependent 
mass More recently, it has been pointed out that the number of 

effective degrees of freedom (i.e. the gluon degeneracy of the system) can 
itself be temperature dependent iHiinj- 

An interesting idea, similar to the one we will use in this work, is imple¬ 
mented in the so-called “cut-off” model cniiiiiiia, in which gluons having 
momenta smaller than a hxed value K are bound inside non-perturbative 
structures and therefore do not directly contribute to the thermodynamics 
of the system. 

In the models discussed so far, two important limitations are present: 
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firstly, the critical temperature plays the role of a parameter and cannot 
be computed within the model; then, and more important, the transition 
itself is parametrized and not obtained dynamically, for instance through the 
minimization of the energy of the system. In the present work we will try to 
overcome such limitations: in our approach it is possible to describe the ther¬ 
modynamical behavior of the system and obtain a hrst order transition via 
the minimization of the thermodynamical potential; the value of the critical 
temperature at which the transition takes place can also be estimated 

In our model, we consider three different contributions to the thermo¬ 
dynamical potential: the hrst component comes from the gluon condensate, 
whose dynamics is expressed in terms of a dilaton lagrangian [TH [T^ [TBl ITT] : 
then, gluons are introduced in a way similar to the one used in the cut-off 
model d Ellis. At variance with the previous versions of that model, in 
our case the cut-off itself is not a parameter, but rather a function of the 
expectation value of the dilaton held, i.e. of the gluon condensate. It seems 
rather natural to assume that, when the gluon condensate is large, namely at 
low temperature, many gluons are frozen inside this non-perturbative struc¬ 
ture: as a consequence, the infrared cut-oh is large. On the contrary, when 
the value of the gluon condensate is reduced, gluons having large enough 
momenta may “evaporate”, and behave as almost-free particles. 

Finally, the last contribution to the thermodynamical potential is due 
to the perturbative gluon-gluon interaction, which is present even at large 
temperature, i.e. for T S> Tc. In the following, we will discuss in sequence 
these three contributions. 


1 Gluon-condensate dynamics 

The idea of a gluon condensate has been introduced many years ago [iHiiini. 
In our approach, we will not try to obtain a mechanism for gluon conden¬ 
sation, namely a model for the QCD vacuum, but we will instead describe 
the dynamics of the gluon condensate by introducing an effective degree of 
freedom, i.e. the dilaton held. 

In fact, there is a deep connection between the gluon condensation phe- 

similar approach has been discussed in Ref. m, where the effective degrees of free¬ 
dom are the dilaton field and constituent gluons which become massive via an interaction 
with the gluon condensate. However, the results obtained in Ref. m are not completely 
satisfactory at temperatures just above the critical one. 
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nomenon and the violation of the scale invariance, which, in QCD at the first 
loop, is qnantihed by: 




( 1 ) 


where Jqcd is the dilatation cnrrent in QCD and is the glnon held 
strength. In order to reprodnce the QCD scale anomaly and to satisfy low- 
energy theorems PI, a dilaton held has been introdnced unmn whose la- 
grangian reads 

V{o), 


Lda = 


where 


D(a) = 


B 


an — a + 4: In 




The violation of the scale invariance is given by: 

r)V 

and it mnst satisfy the equality 

( idil ) — { 3QCD )■ 


( 2 ) 

(3) 

(4) 

(5) 


The potential V{a) has a minimum at a = Uo, where V{ao) = 0; the small 
oscillations around the minimum correspond to the excitations of a scalar 
glueball and they can be parametrized as follows: 


V{a) ~ 2BaQ{a-ao)^ + 0[{a - a^f] 
= \ Mg {a-aof + 0[{a-ao)% 

where the glueball mass Mg is 

Mg = 2 (To ^/B. 


( 6 ) 

(7) 


The dilaton potential contains two parameters, (Tq and S, which can be 
related through eqs. (HHHEUZD to the value of the gluon condensate in the 
vacuum and to the mass of the scalar glueball. Concerning the estimate of 
the gluon condensate (for a review, see 1211 ), Ref. [22] indicates the range 

0.12 GeV^ < {{gGf) < 0.83 GeV\ (8) 
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while both Refs. 0 and M indicate a valne {{gGY) ~ 0.5 GeV"^ associated, 
in Ref. [21], to an error of abont 50%. 

The mass of the scalar glneball has been recently estimated in lattice 
QCD, obtaining a mass Mg = (1730 ± 30 ± 80) MeV ^H] . 

The two parameters cxo and B are therefore constrained into a relatively 
narrow window, the nncertainty being mainly dne to the error bar in the 
estimate of the glnon condensate. 

In the following, we will stndy the thermodynamical potential associated 
with the dilaton lagrangian, at the mean held level, by nsing the standard 
techniqnes of hnite temperatnre held theory 123 EH- 

The thermodynamics of the dilaton held at hnite temperature has already 
been discussed in the literature, for instance in Refs. j2Hl and (221 • 

|22] an attempt to go beyond the mean-held approximation was made, and in 
Ref. [20] the gluon condensate was studied with renormalization group how 
equations. We will compare later our results with the ones of Ref. [201 • 
to the difficulties associated with the quantization of non-polynomial held 
theories (see e.g. ins) we prefer to stick to the mean-held approximation; 
moreover, it is not possible to apply the renormalization group techniques 
to our complete model which incorporates also gluons evaporating from the 
condensate. 

In the mean held approximation the thermodynamical potential reads: 


= V{<t) - Pi„(a,T), 


(9) 


where Pdii, fhe pressure of the dilaton held, reads^ 


Fdil(cr, T) 



In 


e-/T 


( 10 ) 


and 

oj = \Jv‘^ + [m(cr)]2 . (11) 

The cr—dependent mass m{a) is dehned as 



d‘^V 


( 12 ) 


^At very large temperatures the contribution of the dilaton field to the thermodynam¬ 
ical quantities should vanish. To this purpose, in Sec. inm an ultraviolet cut-off will be 
introduced in eq. m- 
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Figure 1: Real part of dilaton thermodynamical potential for various temperatures. 
Temperature increases from upper to lower curves, ranging from 0.2 GeV to 0.4 GeV in 
steps of 0.02 GeV. 


and it equals Mg for a = a^. The mass squared is negative for a < 
and therefore the thermodynamical potential gets an imaginary part for small 
values of a. There is in the literature a broad discussion about the physical 
interpretation of a complex potential and how to deal with it (see e.g. jd2ji. 
We will follow Ref. [SSI, where the imaginary part is interpreted as a sig¬ 
nal of the instability of the system. Our recipe for dealing with a complex 
thermodynamical potential is therefore all simply to minimize its real part. 
Actually the problem of interpreting the imaginary part of the potential, 
although conceptually important, turns out to be not so important from a 
practical viewpoint, since after the introduction of gluons (what will be done 
in the next Section) the dilaton held will vanish for T > Tc and we will not 
have to deal with the instability region. 

In Fig. [U we present Re r2dii(o', T), as a function of a, for various temper¬ 
atures. As one can clearly see, when the temperature increases the real part 
of the thermodynamical potential develops a new minimum for a < (Jq. At 
the critical value T^, the new minimum becomes the absolute one and a hrst 
order transition takes place. We must stress again that the transition of the 
pure dilaton held appears to be hrst order due to the approximations we have 
used. The order of the transition could be established in a consistent way 
only using a non-perturbative approach, e.g. studying the dilaton potential 
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Figure 2: Dilaton field as a function of the temperature in the pure dilatonic model. 

on the lattice. One should also consider the results of Ref. EDI. where no hint 
of a first order transition was found up to temperatures of the order of 200 
MeV where their prediction reaches its limit of validity. On the other hand, 
many calculations of the behavior of the dilaton held at hnite temperature 
do indicate a transition at a temperature similar to the one we get in our 
approach j2Hll2niE21- It is also interesting to notice that, since in the dilaton 
lagrangian the only dimensional parameter is cxo, the critical temperature as 
a function of B and ao must be of the form: 

rf = /(B)<To, (13) 

where f{B) is a function to be determined numerically by minimizing ReOdii- 
In Fig|21 we present the expectation value of the dilaton held as a function 
of the temperature: for T < there is a very small reduction of a from 
its zero-temperature value (Tq. This shift of the expectation value of the 
dilaton held corresponds to thermal excitations of the glueball. At the critical 
temperature the dilaton held is discontinuous. Note that for T > 
the dilaton held does not vanish. We will see later on that the contribution 
of quasi-free gluons shifts the value of a to zero in the deconhned phase. 
Figs. lU and El have been obtained using B = 46.4 and ag = 0.127 GeV, 
which correspond to {{gGY) = 0.35 GeW^ and Mg = 1.73 GeV. The critical 
temperature, in absence of gluons is Tc = = 0.3 GeV: the extra pressure 

of the gluons in the deconhned phase reduces this value, bringing it close to 
the one indicated by lattice calculations. 
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2 The evaporation model 

2.1 Quasi-free gluons 

In the previous section we have discussed the behavior of the dilaton held at 
hnite temperature. It is clear that its excitations (i.e. the glueballs) cannot 
represent the relevant degrees of freedom at large temperature, where quasi- 
free gluons should give the dominant contribution to the thermodynamical 
observables. On the other hand, quasi-free gluons should be suppressed be¬ 
low Tc, which, in our model, is the temperature at which the dilaton held is 
discontinuous. A simple way to suppress the quasi-free gluons in the conhne- 
ment region is by assuming that they are frozen inside the gluon condensate: 
when the value of the gluon condensate is large, i.e. below T^, most of the glu¬ 
ons are frozen while, above T^, the condensate evaporates and gluons become 
quasi-free particles. Technically, this idea can be implemented by introducing 
an infrared cut-oh K in the gluon distribution function, so that only gluons 
having a momentum larger than K contribute to the thermodynamics of the 
system: 


Fq_free(0-,T) = 


-2{Nl - l)T 


d^k 


In 


.-k/T 


Q{k-K{a)) (14) 


In our model we assume that the cut-oh A' is a function of the gluon con¬ 
densate, i.e. of the expectation value of the dilaton held. The “evaporation” 
model has been already discussed in Refs. unmnini, but there the cut-oh 
was assumed to be a hxed parameter. We use for the cut-oh the form 


K{a) 



(15) 


SO that if cr (Tq, then K —> oo, while if a —0 then K ^ A. It can be 
interesting to notice that we can not reproduce the lattice data satisfactorily if 
we use a cut-oh which vanishes for a -C cxo. This result seems to indicate that, 
even at large temperatures, at which a ~ 0, wee gluons are still suppressed. 
The value of the parameter A is of the order of 1 GeV and it is therefore 
natural to interpret this parameter as the one regulating the transition from 
the perturbative behavior of the gluon propagator to the non-perturbative 


one. 







2.1.1 High temperature degrees of freedom 

At very large temperature the perturbative degrees of freedom should be 
recovered and therefore the dilaton held cannot appear as an effective degree 
of freedom for T ^ Tc. On the other hand, at temperatures just above Tc 
scalar gluon-gluon correlations described in terms of the dilaton field could 
still be relevant. We will see in Sec. 3 that, above T^, a ~ 0; as a consequence 
the dilaton mass vanishes and scalar gluon correlations can exist but with a 
vanishing mass gap |2H1- 

The idea of describing correlations between the asymptotic degrees of 
freedom in terms of effective fields is adopted in many physical problems. 
For instance cr-models have been used in studying the chiral phase transi¬ 
tion In that case, the chiral fields describe quark-antiquark correlations 
in an effective non-perturbative way below and above Tc. In the scheme we 
are discussing it is possible, at least in principle, to provide a structure for 
the dilaton field in terms of gluonic degrees of freedom [2]. For simplicity we 
mimic the dynamics of gluons inside the dilaton by introducing an ultraviolet 
cut-off in the dilaton pressure. Since above Tc gluons having a momentum 
larger than K are quasi-free, and since at least two gluons are needed to pro¬ 
duce a scalar correlation, we assume that correlations having a momentum 
larger than 2K are suppressed and do not contribute to the thermodynamical 
quantities. 

We modify therefore eq. (HUD by introducing an ultraviolet cut-off equal 
to 2K. In this way, at very large temperatures, the dilaton field does not 
contribute, as shown in Fig. El On the other hand, the effect of this cut-off 
is almost negligible at T ~ a few Tc. 

2.2 Residual perturbative interaction 

Lattice data clearly indicate that even at a temperature larger than 4 Tc the 
Stefan-Boltzmann limit is not yet reached and the data for pressure, energy 
and entropy he below the free-gas limit. As we will see, the introduction of an 
infrared cut-off is not sufficient to explain both the data near Tc and those 
at large temperatures. It is therefore necessary to introduce perturbative 
corrections, which have a relevant role for T > Tc. These corrections are well 

^Similarly, in the study of hadronic structure, pions can be introduced as effective 
degrees of freedom, whose substructure can later be resolved in terms of the quark distri¬ 
bution function of the pion (see for instance Esiisni)- 
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established and a huge amount of work has been devoted to their estimate 
(for a recent review see [SZ])- In fhe present work we are not interested 
in a detailed comparison with the data at T S> Tc, but we mainly aim at 
describing the data near the critical temperature. We have therefore adopted 
the simplest prescription, which consists in taking into account only the hrst 
order 0{g‘^) corrections. These corrections have been computed in Ref. jd8j . 
but in the present case additional ©-functions occur which restrict the phase 
space for the interacting gluons (see 0). The 0((7^) contribution to the 
pressure reads: 


Pi,t(a,T)=/W(iV'-l)|-3(^| 

f d^ki 


d^/c2 1 

(27r)^ J (27r)^ kik2 

'9 


ki 


d^k 1 
(27r)3 k 
'h 


N, 


B 


T 


Q[k-K{a)] 


iVs ( ^ ) iVs ( ^ ) Q[k, - K{a)] Q[k2 - K{a)] 


T 


X 


0[|ki + ksl - K{a)] - - 0[|ki - ksl - K{a)] 


(16) 


where iVs(a:) = (e^ — 1) ^ is the Bose-Einstein distribution. Here is the 
temperature-dependent running coupling-constant 


9\T) 


487r2 

lliVdn[(T2 + p2)/A2] ’ 


(17) 


where we have introduced a regulator S', whose appearance can be related to 
the existence of a minimal momentum K for the propagating gluons [3]^. As 
we shall see, the perturbative corrections play a relatively minor role in our 
model and more sophisticated choices of the running coupling would hardly 
affect the results. The typical value for S is S ~ GeV. 


principle, g'^ can also depend on the value of the gluon condensate PHI- This 
possibility will be explored in a future work HI!- 


10 








3 Results 

The parameters in our model are the following: 

• dilaton lagrangian: B, cxo 

• quasi-free gluons: A, a 

• running coupling-constant: A, S. 

Concerning B and (Tq, as discussed in Sec. 1 their value is bounded by the 
“experimental” value of the gluon condensate and by the lattice result for 
the mass of the scalar glueball Mg. In the following we will present results 
obtained using B = 46.4 and do = 0.127 GeV, which correspond to {{gG)"^) = 
0.35 GeV^ and to Mg = 1.73 GeV, both near the center values indicated for 
these quantities. Goncerning the quasi-free gluons, we use A = 1.015 GeV 
and a = 0.5. Finally, the parameters for the running coupling-constant 
are not too strictly constrained in our calculation, since the perturbative 
interaction turns out to play a minor role in our model. We adopted S' = 7.15 
GeV and A = T^, where is the value of the critical temperature computed 
in the model. These parameter values are also consistent with the more 
general form for introduced in HOI (see note ID). Notice that the value of 
the critical temperature is not modihed by the perturbative corrections and 
can therefore be computed before the latter are taken into account. 

3.1 Thermodynamical quantities 

The results we present are obtained by minimizing the real part of the total 
thermodynamical potential Gtot(c’‘, 1"), as a function of a, for a given temper¬ 
ature T, where: 

Gtot(o-, T) = Gdii(cr, T) - Pq_free(cr, T) - Pint(cr, T) . (18) 

In Fig. 121 we present the real part of the thermodynamical potential. Due 
to the extra pressure of quasi-free gluons, a new minimum develops for cr = 0. 
In this way we avoid entering the region of instability of the dilaton held, 
where imaginary parts develop, since the expectation value of the held jumps 
from a value a ~ cxo to a = 0. The value of the critical temperature, which 
corresponds to the discontinuity of the dilaton held, is = 0.27 GeV. The 
reduction of the value of the critical temperature, due to the introduction 
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a [GeV] 


Figure 3: Real part of total thermodynamical potential for various temperatures. Tem¬ 
perature increases from upper to lower curves, ranging from 0.23 GeV to 0.33 GeV in steps 
of 0.02 GeV. 

of gluons in the deconfined phase, is rather independent of the precise value 
of the cut-off parameters and it turns out to be always of the order of 10%. 
The critical temperature is therefore a rather good approximation to the 
value of Tc as computed from the complete model. 

In Fig. m we compare our result for the scaled pressure with the lat¬ 
tice data PJ. The pressure is obviously connected to the thermodynamical 
potential by the relation P = —hltot- We also show the pressure obtained 
minimizing the thermodynamical potential neglecting the interaction term 
Pint- It is clear that the interaction modifies significantly the pressure only 
for temperatures of the order of 2Tc, or larger. The critical temperature, 
as well as the shape of the pressure near T^, are independent of the inter¬ 
action contribution. As already stated, we are not particularly interested in 
reproducing in a very accurate way the data at large temperatures, introduc¬ 
ing perturbative corrections in a sophisticated way, our main aim being to 
describe the phase transition. 

In Fig. 13 we present the decomposition of the total pressure into its vari¬ 
ous contributions. As it can be seen, the dilaton contribution and the quasi- 
free-gluons one are both discontinuous at the critical temperature, but their 
discontinuities cancel so that the total pressure is continuous. In Ref. [2], 
such a behavior for the gluon condensate and for the quasi-free gluons was 
anticipated. From the figure it is also clear that the perturbative correction 
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T/Tc 


Figure 4: Scaled pressure in our model compared with lattice data. The solid line is 
obtained minimizing the total thermodynamical potential, while the dashed line is obtained 
neglecting Pi„t. 
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Figure 5: Various contributions to the total pressure. The long-dashed line corresponds to 
Pdii = — Udii, the short-dashed line is Pq_free, the dashed-dotted line is Pint and the dotted 
line is Pdii but without the ultraviolet cut-off. 
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Figure 6: Scaled energy density in our model compared with lattice data. The solid 
line is obtained minimizing the total thermodynamical potential, while the dashed one is 
obtained neglecting Pint- The shaded rectangle indicates the latent heat as estimated in 
lattice calculations. 

Pint vanishes at T^. Let us remark again that the dilaton gives a contribu¬ 
tion to the scaled pressure which vanishes at large temperatures, due to the 
presence of the ultraviolet cut-off. 

In Fig. El we compare the energy density e = — P computed in our 

model with lattice data®. We also show the lattice result for the latent heat. 
The main difference between our result and the lattice one is that at T < Tc 
our energy density is considerably smaller than the one indicated by lattice 
calculations. This discrepancy is due to the presence, in our calculation, of 
the scalar glueball only, while the J = 2 glueball should also contribute. The 
mass of the latter has been estimated to be M 2 ++ = 2400 ±25 ±120 (221 • The 
introduction of these new degrees of freedom would correspond, roughly, to 
a degeneracy factor of 6 in front of the glueball contribution [221; and would 
bring the computed energy density close to the lattice one at T < Tf.. For 
simplicity we have not included the excitations of the tensor glueball in our 
calculation, but it can obviously be done in the future. 

We must also notice that lattice data for the energy decrease much faster 
than the results of our model for temperatures just above T^. We have found 

®We computed also the entropy density in our model. Considerations similar to the 
ones done for the energy density can be done for the entropy. 
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Figure 7: Interaction measure. The solid line is the total value, the dot-dashed is the dila- 
ton contribution, the dotted is the one due to quasi-free gluons and the dashed corresponds 
to the perturbative interaction. 


that, using other parameter sets, it is possible in principle to reduce this 
discrepancy. However, in those cases the dilaton field does not jump directly 
to zero at Tc, but it reaches a small but finite value. Due to the difficulties 
associated with the mean field treatment of the model for small values of a 
we have not explored this possibility in detail. 

Finally, in Fig. [7|we show the various contributions to the “interaction 
measure”, namely the (scaled) quantity indicating the distance from the 
Stefan-Boltzmann relation e = 3P. For T T^. the main contribution to 
this “interaction measure” comes from the dilaton field. The contribution 
due the quasi-free gluons is large at moderate temperature, but it decreases 
more rapidly than the contribution due to Pint, which is the dominant term 
at very high temperature. 
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3.2 Thermal gluon masses 


In a covariant gauge, the gluon propagator can be written in the following 
general form 


D^u = 



pr + 


1 







(19) 


where P^^ and P^^ are the longitudinal and transverse projection operators, 
defined as 


po* = Pi? = 0 

- P^^ (20) 

The gluon self-energy is given by 

= GP^'" + FPl\ ( 21 ) 

G and P being scalar functions of q^ and |q|. The electric and magnetic 
masses are dehned as 


pOO 





P(0,q^0) = -noo(0, q-^ 0) =771^; 

G(0,q^0) = ^n,,(0,q^0)=mLg, (22) 

where the relation between P, G and 11 comes from eqs. (EOl) and (ED). 

The gluon self-energy can be evaluated perturbatively; we will concentrate 
on the electric mass, due to the difficulties associated with the magnetic one. 
Two different contributions arise for n°°, a zero temperature one, which 
vanishes after renormalization, and a finite temperature one, which reads 


noo 



d^p 1 1 

(27r)^ uj — 1 


(23) 


In our model, which allows contributions to the thermodynamics of the sys¬ 
tem only from gluons having a momentum larger than K, we rewrite the 
above formula as: 


noo 



d^p 1 1 

(27r)^ UJ e^^^ — 




(24) 
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Figure 8: Gluonic electric mass in our model (solid line), and in a purely perturbative 
calculation without cut-off (dashed line). The data, taken from Ref. [IT], correspond to half 
the screening mass /r(T)/2 (see text). The dotted line indicates niei in the metastable phase. 


In Fig. IHl we show our results; the gluonic electric mass presents a disconti¬ 
nuity at T = Tc, because of the discontinuity of the a held. We also show 
(dotted line) rriei in the supercooled, metastable phase in which a is kept 
equal to a(T = T+) for T <Tc (see page. 152 of Ref. |2n|)- For comparison 
we present the electric mass in a standard perturbative calculation without 
the cut-off, which yields mei = gT (dashed line). 

To compare our results with lattice QCD calculations we must recall that 
a screening mass fi{T) can also be introduced by considering the behavior 
of the potential V{R, T) between gauge invariant sources at T > T^. This 
potential can be parametrized as [¥!] : 


V{R, T) _ e{T) 
T ~ {RTY 


(25) 


where d, e(T) and /i(T) are determined from lattice results. In perturbative 
calculations, the screening mass p-(T) and the gluon electric mass niei are 
connected by the simple relation 


yu = 2mei . 


(26) 
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In Fig. IHl we see that this perturbative relation is recovered for T ^ T^. For 
temperatures just above T^. the situation is less clear, also due to the large 
error bars in lattice data, but the simple proportionality indicated by eq. (ESD 
seems to be (slightly) violated. 

Concerning the magnetic mass, it is well known that it vanishes in a per¬ 
turbative calculation. Notice anyway that, due to the infrared cutoff which 
characterizes the evaporation model, higher order loops are not expected to 
be affected by infrared divergences, even in the absence of a magnetic mass. 
For the same reason, in a cut-off model like the one we have discussed here, 
perturbative corrections can be computed in principle up to an arbitrary 
order; moreover, all these contributions vanish at T = T+. 


4 Conclusions 

We have presented a simple realization of an evaporation model, in which 
at low temperature gluons are frozen inside the non-perturbative condensate 
while at high temperature they escape from the condensate and behave as 
quasi-free particles. We have shown that it is possible within the model to 
reproduce the main results of lattice QCD for thermodynamical quantities 
such as pressure and energy. At variance with other models for these quan¬ 
tities, in our approach a hrst order transition is obtained by minimizing the 
thermodynamical potential and the latent heat can be estimated. It is also 
possible within the model to study hnite temperature gluon masses. Here 
again our results are consistent with the indications of lattice QCD. 

The main problem in the present analysis arises from the difficulties as¬ 
sociated with the appearance of an imaginary part in the thermodynamical 
potential. We assumed that the imaginary part signals an instability of the 
system and we have therefore minimized the real part of the thermodynam¬ 
ical potential. This difficulty is clearly related to the mean-held approxima¬ 
tion we have adopted. Although this technique seems sufficient to reproduce 
rather precisely lattice QCD results, it is clear that only more sophisticated 
approximations can clarify the details of the behavior of the system near the 
critical temperature. However, in the complete model, which includes gluons 
in the deconhned phase, the expectation value of the dilaton held is such 
that the system does not enter the unstable region at any temperature. 

This calculation can be extended in several directions. Firstly, it will 
be important and, probably, relatively easy to study the behavior of the 
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phase transition as a fnnction of the color nnmber N^. Lattice calcnlations 
discnssing the dependence on Nc of the critical temperatnre, of the glne- 
ball masses and of the glnon condensate appeared recently (see for instance 
Refs. j42114dj i . allowing comparisons with onr model. As pointed ont in (2], 
in this model the valne of the critical temperatnre is Nc independent in the 
large Nc limit. 

Another, more important extension of the present work will be the in- 
clnsion of quarks. This can be done by dressing the quark propagator via 
the Schwinger-Dyson equation and by using a quark-gluon coupling which 
depends on the value of the gluon condensate. Work along these lines is now 
in progress m- 
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M.P. Lombardo, R. Tripiccione for many useful discussions and F. Karsch 
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